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The decay of metastable patterns in the form of magnetic hquid ridges is studied in the 
frame of a hnear stability analysis where the ridges are the most unstable linear pattern of 
the Rosensweig instability. The analysis of the corresponding dispersion relation reveals 
that different sets of solutions exist. Two of them are associated with an oscillatory and 
a pure exponential decay, respectively. 

1. Introduction The most well known phenomenon of pattern formation 
in magnetic fluids is the Rosensweig or normal field instability. This instability can 
be observed for a horizontal layer of magnetic fluid (MF) with a free surface with air 
above. Beyond a certain threshold Be of a vertically applied homogenous magnetic 
induction, the surface becomes unstable and a stable pattern of peaks develops [1]. 
That final pattern, resulting from nonlinear interactions, is one possible state 
towards which the preceding most unstable linear pattern of liquid ridges can 
develop if a supercritical induction Bsu-p > is abided [2]. The other possible 
state results from the decay of the liquid ridges towards the flat surface if the 
induction is switched from a supercritical to a subcritical value i?sub [3]. With 
these two possible final states, the most unstable linear pattern can be considered 
as a metastable pattern formed by magnetic liquid ridges. 

That ridges were studied in [3,4] by means of an experimental set-up in which 
the magnetic induction was jumplike increased from a start value of < B^ 
to i?sup [4] and after a short time delay jumplike decreased from Bsup to Bsub 
[3]. The jumplike change of the magnetic induction allows to apply the linear 
stability theory, where the magnetic induction is assumed to be instantly present. 
Consequently, the quantitative agreement between measured and calculated data 
is very good for the wave number of the linearly most unstable ridges [4] and for 
the propagation velocity as well as the oscillation frequency during the decay of 
the ridges [3]. The calculation of the latter two quantities disclosed new features 
of the dispersion relation which will be revisited and analysed in the next section. 

2. Analysis of the dispersion relation A horizontally unbounded layer 
of an incompressible, nonconducting, and viscous MF of infinite thickness and 
constant density p is considered. The MF has a free surface with air above, where 
the basic state is that of a nondeformed surface. The dispersion relation for small 
disturbances from the basic state is given by [5] 
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where g is the gravity acceleration, u the kinematic viscosity of the MF, a its 
surface tension with air, p,r its relative permeability, B the absolute value of the 



external magnetic induction, and /io the permeability of vacuum. For the dis- 
persion relation (1), the small disturbances were decomposed into normal modes, 
i.e. they are proportional to exp[— — qr)]. r = {x,y) is the planar space vec- 
tor, q = (qx,qy) the wave vector, and q — |q| denotes the wave number. With 
u) = uii+iLV2 the real part of —iiv, UJ2, is called the growth rate and defines whether 
the disturbances will grow {u>2 > 0) or decay {u2 < 0). The absolute value of 
the imaginary part of — iw, \uJi\, gives the angular frequency of the oscillation if 
Ui is different from zero. The critical induction for the onset of the Rosensweig 
instability at w = is = 2/xo/Xr-(A<r + 'i-)y/p(Tg{^r — [!]• 

Figure 1 shows for the reasons of clarity one solution of the dispersion relation 
for a supercritical induction Bgup = 1.05i?c (thick solid line) and for a subcriti- 
cal induction Sgub = 0.96Sc (thick dashed line). In the supercritical case one 
can conclude from Re[uj^{qyn, Bsup)] = ^i.m ^ ^l.m < Ini[Lt^^((Zm, -Bgup)] = 
2a; i^m 1^2,111 = that a most unstable linear patterns evolves with the maximal 
growth rate uim = *W2,m (w2,m > 0) and its corresponding wave number [thin 
solid lines in Fig. 1(a)]. In the subcritical case the relations Rc[a;^((7ni, -Bsub)] > 
and Im[w^(gm, -Bgub)] < [thin dashed lines in Fig. l(a)(b)] do not allow a simple 
conclusion on the values of wi and i02- Therefore the classical presentation of 
versus q chosen for the analysis of inviscid MFs [6-8] is only of limited use for the 
analysis of viscous MFs. 

A more meaningful display is the plot of Re(— iw) and Im(— iw) versus the 
subcritical induction for different viscosities (see Fig. 2), where (jm(i^) of the most 
unstable linear pattern was determined for Bg^p = 18 mT. If the fluid is inviscid 
(dashed lines in Fig. 2), the pattern oscillates below a certain threshold for the 
subcritical induction, -Bc.o['?m(i^ = 0),cj2 = 0] [sec Eq. (5)]. Above that threshold 
the pattern can either decay according to the solution Re(— iw) < or can develop 
towards the most unstable linear pattern belonging to Bg^^ > Be since Re(— iw) > 
is also a solution. 

For a viscous fluid (solid lines in Fig. 2) the behaviour is more complex. A 
first critical induction, Bc,\, occurs, where the set of solutions for the dispersion 
relation (1) changes from two complex solutions to two negative real solutions. 
[The same transition occurs for the damping of gravity-capillary waves as param- 
eters as the aspect ratio and the so called inverse Reynolds number are varied 
(see Figs. 4 6 and 8,9 in [9]).] Both real solutions exist imtil at a second critical 
induction, iJc,2, one of them abruptly ends. At a third critical induction, B^fi, one 
of the two negative real solutions changes its sign and becomes positive. 

To understand this complex behaviour, the dispersion relation (1) in dimen- 
sionless units (indicated by a bar) is analysed in the rearranged form 

-iu!\^ ^ -iw -q - q^ + 2B'^q^ 

• (2) 
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All lengths were scaled with [(7/(^(7)]-'^/^, the time with a'^/^ / {g'^/^ p^/'^), the vis- 
cosity with cr^/^/((/^/^p^/'^), and the induction with Be- Equation (2) reveals that 
whatever value (D G C has, the left hand side of Eq. (2) has to be real because 
the right hand side of Eq. (2) is always real since (g, P, B) G M. That condition 
together with the mixing of real and complex quantities in Eq. (2) is essential to 
understand the above described appearance of different sets of solutions. As long 
as Bsub ^ Bc.i all solutions of the dispersion relation with q = qm are real, i.e. 
—iw = a;2 e K (see Fig. 2). Using this result it follows from Eq. (2) that there is a 
value u)2 = —i^Qm beyond which the radicand becomes negative. Since a complex 
value for the left hand side of Eq. (2) is not allowed, the solution does not exist 
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beyond uJ2 = ^'^Qm- This corresponds to the point in Fig. 2, where one of the so- 
lutions suddenly terminates at i?c,2- Therefore the second critical induction yields 
to 
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The first critical induction is the minimal induction for which real solutions exist, 
thus 



-Bc,l[^>9m(l^),W2,r 



\ 2 U 



VTTd- ( 1 + 



D 



(4) 



witha}2,min = Pql,D = [- 60+ (108+ 12^93)2/3 + 144(108+ 12 V93)-2/3] . 

Finally, the third critical induction is defined by OJ2 = which leads to 



Bc,o[qin{i^),i02 = 0] = y ^ i^Y' 



(5) 



At a closer inspection of Eqs. (3-5) one realizes that the three thresholds follow the 
relation Bc,i < Bcfi < -Bc,2- Their dependence on the viscosity of the magnetic 
fluid will be studied in the next section. 



3. Results and discussion The behaviour of the three above defined 
thresholds on the viscosity is studied, where all other material parameters are 
kept constant. The wave number gm(z^) for the most unstable linear pattern is 
determined for a supercritical induction of 18 mT. Figure 3 shows that for viscous 
MFs four regions with different sets of solutions for the dispersion relation exist. 
For small subcritical inductions the set of solutions consists of two complex solu- 
tions with negative real parts (symbolised by 2C__ in Fig. 3) which is followed 

by two negative real solutions (2R ). Passing the threshold i?c,o (solid line), the 

set of solutions is formed by one positive and one negative real solution (2R_| ). 

Only the former solution lasts at high subcritical inductions (1R+). From Fig. 3 it 
becomes clear that for low viscous MFs, D < 0.1, only an oscillatory decay can be 
observed [3, 10]. For high viscous MFs it should be possible to observe an oscilla- 
tory (2C__) as well as a pure exponential decay (2R__) of the pattern. Because 
for > 0.3 the two thresholds Bd (dashed line) and i?c,o are well separated and 
the region 2R becomes experimentally accessible. 

For a nonmagnetic fluid with the same density and surface tension, a dimen- 
sionless viscosity of i> ~ 0.9 is necessary to have a pure exponential decay, i.e. the 
nonmagnetic fluid has to be considerable more viscous than the magnetic fluid to 
result in the same behaviour. Another advantage in the use of a MF is that the 
jumplike increase to -Bgup is a simple way to prepare just one particular metastable 
pattern whose decay can be studied. 

For inviscid MFs Fig. 3 shows that only two different sets of solutions occur. 
For subcritical inductions smaller than Bcfi{y = 0) two pure imaginary solutions 

(21) exist. Above that threshold two real solutions with different signs (2R_| ) are 

found. 

This analysis also leads to a better understanding in the case of a supercritical 
induction, where regions with different sets of solutions should exist, too. Whereas 
Fig. 1 shows only one particular solution of the dispersion relation (1) for < 
q < 1200 m-\ the two zooms around 5=418.3 m'^ [Fig. 4(a) (b)] and g=1027 
m~^ [Fig. 4(c) (d)] present all solutions. In contrast to the material parameters 
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for Fig. 2, now the induction is fixed at -Bsup = 18 mT and therefore transitions 
between different sets of solutions occur at certain wave numbers. The first zoom 
[Fig. 4(a) (b)] displays that for wave numbers below gc,i two complex solutions 
exist. Between cjci and gc,2 two real solutions exist and beyond (jc.2 only one 
real solution. In a reverse order the different sets of solutions are appearing in 
the second zoom [Fig. 4(c) (d)]. These features of the dispersion relation have not 
been realized since only the solution for was of interest in [2,4, 11]. For that 
particular wave number the dispersion relation has only one real solution which 
motivates the plot of one solution (thick solid line) in Fig. 1. 

To conclude, in the frame of a linear stability analysis the decay of the 
metastable pattern of magnetic liquid ridges is studied. There are initiated by 
a jumplike increase of the magnetic induction to a supercritical value whereas the 
decay is triggered by a jumplikc decrease of the induction. The analysis of the 
dispersion relation according to such a procedure reveals that depending on the 
value of the subcritical induction different sets of solutions exist (Fig. 2). The 
regions for these different sets are separated by the three thresholds (3-5) whose 
dependence on the viscosity is displayed in Fig. 3. Whereas for low viscous MFs 
an oscillatory decay of the ridges can be observed [3] , for high viscous MFs a pure 
exponential decay should be detectable, too. 
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Figure 1: Dependence of Re(a;^) (a) and Im(a;^) (b) on the wave number q for 
a supercritical (1.05i?c, tliick solid line) and a subcritical {0.96Bc, thick dashed 
hne) induction. The wave number of the most unstable linear pattern and 
the corresponding value are indicated by thin solid lines. The solution in the 
subcritical case for q^ is indicated by thin dashed lines. The material parameters 
for the calculations are u = 5.17- IQ-^ mVs, p = 1.16- 10^ kg/m^, a = 2.65-10-2 
kg/s^, iir ^ 1-935, and Be = 16.84 mT. 
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Figure 2: Dependence of Re(— iw) (a) and Im(— iw) (b) on the subcritical induction 
Bsub for an inviscid MF {u = 0, dashed hne) and a viscous MF {v = 5.17 • 10~^ 
m^/s, solid Hnc). For the former fluid two imaginary solutions arc replaced by two 
real solutions with increasing -Bgub- For the latter fluid a region with two complex 
solutions is followed by a region with two real solutions which is succeeded by a 
region with one real solution. The three critical inductions -Bc,0! Sc,i: ^^nd _Bc,2 
are explained in the text. The remaining material parameters for the calculations 
are those from Fig. 1 with qni{i^) determined for Bgup = 18 mT. 
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Figure 3: The three critical inductions Bc,i (dashed line), Bc,o (solid line), and 
Bc^2 (long dashed line) versus the dimcnsionlcss viscosity D. The regions with 
different sets of solutions for the dispersion relation are symbolised as follows: 21 - 

two imaginary solutions, 2C - two complex solutions with negative real parts, 

2R_| - two real solutions of different signs, 2R - two negative real solutions, 

1R+ - one positive real solution. 
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Figure 4: All solutions of Eq. (1) for B = 1.05Sc in the form of Rc(ti)^) (a)(c) and 
Ini(ti;^) (b)(d) for wave numbers near q— 418.3 (a)(b) and 9=1027 m^^ (c;)(d)- 
In (a)(b) two complex solutions for q < ^c,! are followed by two real solutions 
which exist for Qc^i < q < 9c, 2 (dashed and long dashed lines). For (?c,2 < q one 
real solution (long dashed line) exist. In (c)(d) the order of the different solutions 
is reverse: first one real solution, then two real solutions and finally two complex 
solutions. The material parameters are those from Fig. 1. Note the different scales 
at the axes. 
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